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REVISED ABSTRACT 

Many calculations of the electrical conductivity of nonequilibrium 
MHD plasmas have been based on a simple two- temperature theory. This theory 
assumes that (a) the electron number density has the equilibrium (Saha) 
value corresponding to the electron temperature T e , and (b) the free elec- 
tron energy distribution function f(u) is maxwellian. 

The validity of assumption (a) has been studied by several authors 
who, however, used assumption (b) in their analyses. Assumption (b) has 
been less thoroughly studied. Because both excitation and ionization rates 
are sensitive to f(u) and bound states may be out of equilibrium due to 
radiative transitions, it is unrealistic to treat these two assumptions 
separately. This paper reports preliminary results of an investigation 
undertaken to establish the range of validity of the two-temperature theory 
for MHD plasmas by solving the Boltzmann equation for f(u) and the steady- 
state rate equations for the bound electronic states. The problem was 
attacked in the following three stages. 

First, f(u) was calculated from a Boltzmann equation including only 
the electric field and the elastic collision terms. The results showed 
that for typical MHD systems (e.g., Ar + K at 1 atm) electron-electron 
collisions drive f(u) to maxwe Ilian. 

Second, the solution of the rate equations for a maxwellian f(u), 
using a five-level cesium atomic model, demonstrated the importance of 
radiative transitions in determining the bound-state populations and 
magnitudes of the inelastic collision terms. The model atom consisted of 
four discrete states; (6S; 6P,P'; 5D,D'; 7S) and a lumped state, which was 
assigned various binding energies and degeneracies. Criteria for se- 
lecting the latter were based on the maximum stable orbit radius that 
would be likely for the plasmas of interest. Both the classical Bohr- 


Thomson and Gryzinski cross sections were used to calculate the rate coef- 
ficients and collision terms for excitation, deexcitation, ionization, and 
three-body capture. Many calculations were also performed using recent 
experimental data for the 6S -» 6P,P* excitation. Known values of radia- 
tive probabilities were used for transitions between discrete states. 

Finally, the simultaneous solution of the Boltzmann equation (including 
all relevant elastic and inelastic collision terms) and the bound-state rate 
equations was attempted by an iterative technique. Plasmas of various 
optical thicknesses to resonance radiation were considered. In the extreme 
case of a gas transparent to all line radiation, the departure of the state 
populations from equilibrium causes a large reduction in the free electron 
number density, especially at low values of T e . Because excitation colli- 
sions are not balanced by superelastic collisions in this case, they cause 
f(u) to depart significantly from maxwellian at high energies. Numerical 
difficulty was encountered in the iterative scheme in this case. The 
effect of inelastic collisions is much less pronounced for plasmas that 
are optically thick to resonance radiation, and it disappears altogether 
for a plasma optically thick to all radiation. The distribution function 
was calculated for a range of electric fields and gas densities for the 
case where resonance radiation is completely trapped, the most realistic 
situation for large MHD generators. The dependence of the electrical con- 
ductivity on the current density was also calculated, and where possible 
compared with experimental results. 
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INTRODUCTION 

A fundamental problem of gas discharge theory is the prediction of the 
free electron number density and velocity distribution function, which 
together determine the electrical properties of the gas. In the case of 
certain low-voltage atmospheric -pres sure glow discharges in alkali-seeded 
noble gases, which are being studied for nonequilibrium magnetohydro- 
dynamic (MHD) power generation, this problem is fortunately somewhat 
alleviated. For such discharges, a simple two -temperature conduction 
theory due to Kerrebrock [l] has shown reasonable agreement with experi- 
ment at sufficiently high current densities (above 10 4 amps-m -2 ). 
Kerrebrock 's theory assumes that (a) the free electron energy distribu- 
tion f(u) is maxwellian, with a temperature T e that may considerably 
exceed the gas temperature T g , and (b) the free electron number density 
N e is given by the Saha equation applied at temperature T e . 

There remains, however, the question of the range of validity of the 
two- temperature theory. Kerrebrock and Hoffman [2] attributed the lack 
of agreement between theory and experiment at low current densities 
(below 3X10 3 amps-m -2 ) to the fact that electron-electron collisions 
could not preserve a maxwellian velocity distribution in the presence of 
elastic electron-atom collisions, due to the low value of N e . They 
supported their hypothesis by estimating average frequencies of energy 
transfer through these collisions. 

Archambault and Cuenca [3], on the other hand, concluded that the 
distribution function would be maxwellian under MHD conditions. This con- 
clusion was based on a detailed comparison of the various elastic and in- 
elastic collision frequencies in an argon-cesium plasma, unfortunately for 
only one set of values of N e , T e , Tg, and argon and cesium number densi- 
ties. The authors of [3] raised doubts concerning the validity of the 
Saha equation, however - a position which is contrary to the conclusions 


of an earlier study by Ben Daniel and Tamor [ 4] . 

The present work attempts to ascertain the validity of the two- 
temperature theory for alkali-seeded argon MHD plasmas by calculating f(u) 
and N e self-consistently from the Boltzmann and rate equations. Such a 
self-consistent calculation is necessitated by the fact that f(u) depends 
strongly on N e [5], which in turn is very sensitive to the behavior of 
f(u) above the excitation and ionization thresholds [6]. 


THEORY 

Electron Boltzmann equation 


• The energy distribution function-*- f(u) is defined such that 
u l/2f(u)du is the fraction of electrons having energies in the range from 
u to u + du electron volts; hence, 


f u 1 /2f( u )du = 1 (1) 

J 0 


For a steady-state discharge in a uniform, partially ionized gas flow- 
ing through crossed electric and magnetic fields, the appropriate equation 
for f (u ) is 2 [5,7,8] 


(E*) 2 Mu) 

3 1 + CD 2 * 2 



h=a,i 


+ 


kT g df 
e du 


(i) 


(ii) 



N e u 2 Q“ e (u) 



( u * )3/2 f ( u . ) du . 


+ 



f ( u 1 )du 


— * 

t 


( iii) 


+ f(u) f (u' ) 1 / 2 f(u' )duj = I(f ) (2) 

(iv) 

Equation (2) is sometimes called the gain equation, because when 
multiplied by N e (2e/m e )l/2 its terms give the rates (in m - 3-sec"l) 
at which electrons are added to (or lost from) the energy range u to oo. 
The terms are as follows : 

(i) Energy gain from the effective electric field US* 


^Strictly speaking, u-*-/ 2 f (u) is the energy distribution function, 
but the above terminology will be used for brevity. 

2 Symbols are defined in appendix A. 
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(ii) Energy loss to atoms and ions in elastic collisions 

(iii) Energy exchange in electron-electron collisions 

(iv) Energy gain (loss) in inelastic (including superelastic ) 
collisions with atoms and ions 

The inelastic gain term 1(f) is the sum of the following expressions 
obtained by integrating the inelastic collision terms given in [8,9] from 
u to oo, making simple transformations, and using the microscopic 
reversibility relation between excitation and deexcitation cross sections 
(Klein-Rosseland relation) and between ionization and three-body capture 
cross sections [9]: 

(a) Excitation and deexcitation by electron impact 



X u-cg K (u-)au' (3) 


(b) Ionization by electron impact 


I ion^ f ^ 




u'f(u' )du' 


/ u'-E L 


4 on (u",u') 


*1 


+ af on (u' - u" - E l ,u' ) 


du" 


/•“ /• U '" E I 

- I u'f(u')du' / 

•(i 


al on ( u",u')du"> (4) 


( c ) Three-body ion-electron-electron capture'' 


E cap( 


= h2 , ) 

2 \2«m e e ) 


3/2 


N e N 




u'-E l 


f(u")f(u' - u"- E l ) 


X a3' on (u ,, ,u' )du" 


>•00 /•u'-Ej 

- I u'du' / 

*Ve t •'u 


f ( u" ) f ( u ' - u" - E l ) 


+ a£ on (u’ - u" - E l ,u 


p^ on (u",u' 

'T't 


(5) 


'-’Three-body collisions in which the third body is an atom or ion will 
not be considered. 
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(d) Radiative capture 4 


u'f(u' )o£" ca P(u' )du' (6) 

L 

When u = 0, the excitation-deexcitation term (3) vanishes because 
Q£^ k ( u ) = 0 for u < £E l k , while the terms (4), (5), and (6), when 

multiplied by Ng^e/mg) 1 / 2 , become the total ionization and capture rates. 
Also, when the bound and free electrons are in equilibrium, the term (3) 
vanishes while (4) and (5) cancel. 



Rate equations for bound electrons 

The rate of change of the number of bound electrons in level L is 


N t = -N t (N 3{ + A) + N 

L h e L L e 


K?tL 


Vk-L + I Vk + + P L ) 


(7) 


K>L 


For the steady-state problem treated herein, N L = 0. Equation (7) together 
with the normalization condition 






( 8 ) 


and the condition of charge neutrality (N e = Nj.) serve to determine N e 
and Nl for a specified f(u) and initial seed number density N^. 

Electrical conductivity 

The electron current density j in the plane perpendicular to B 
is [7] 

j = CjE* + a H | X E* 

where 


( 9 ) 


1 « /2e\ 

a. = ■=■ N e — 

1 3 e / 


®f /' rrte “">"(■ 8- 


( 10 ) 


°H = 


i H e e(2l') l/2 r ^ - MuW- SlW 

3 e \ m e / J l + a) 2 x 2 V du; 


( 11 ) 


are the perpendicular and Hall conductivities, respectively. 


^The inverse process to (d), photoionization, is ignored. 
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Electron energy equation 


The electron energy balance equatioft may be derived from equation ( 2 ) 
by multiplying by K e e(2e/m e ) 1 / 2 and integrating from u = 0 to oo. 
After considerable transformation, the result may be written as 


j 


E* 



+ ]^L,K(Vlr>K ~ %%*) + 2 E L( N I^L 0n - Vh aV ) 

L,K>L L 



( 12 ) 


Equation (12) can, of course, also be derived from elementary considera- 
tions. By multiplying equation (7) by E L and summing over all bound 
levels L, one can show that the inelastic terms in (12) (the last three 
on the right) are equal to the power loss per unit volume through line and 
continuum radiation; that is, 



Z I 

L,K<L 


e ^ K ,L N L A L-*K + N e N i e 



(E^ + u)uf (u)o£" Ga P(u)du^iLj 


(13) 


For a maxwellian f(u), the elastic loss term (the first on the right 
of (12)) is proportional to the temperature difference T e - T g , and equa- 
tion (12) provides a means of calculating T e . If f(u) is no? maxwellian, 
the temperature can be determined from the relation 



(14) 


CALCULATIONS AND RESULTS 
Computational procedure 

An iterative technique was employed to solve equations ( 2 ) and ( 7 ) . 

The coefficients of the rate equations (7) are first calculated for an 
assumed f(u) (usually a maxwellian with a guessed T e ). Solution of these 
equations yields N e and Np. When these number densities and the 
assumed f(u) are used, the inelastic term 1(f) and the integrals in the 
electron-electron collision term are then calculable. This reduces the 
equation for the next approximation to a first-order linear differential 
equation in f(u), which is integrated numerically. The solution is then 
normalized to satisfy equation (l) and an electron temperature calculated 
from equation (14). When f(u) is very close to maxwellian, equation (14) 
becomes an identity, and it is necessary to determine T e by satisfying 
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the energy equation (12). Otherwise, the process is repeated until con- 
vergence is obtained in f(u), T e , and N e . 


Results of elastic -collision problem 


Initial calculations were concerned with determining f(u) from equa- 
tion (2), neglecting the inelastic term 1(f). These calculations were 
performed for argon seeded with potassium and cesium under the conditions 
of Kerrebrock and Hoffman’s experiment, to see whether a non -maxwe Ilian 
distribution would be obtained at low current densities, as hypothesized 
in [2]. For these calculations there was no magnetic field, and N e was 
assumed to be given by the Saha equation. Elastic cross sections for 
argon, potassium, and cesium were taken from [10, 11, 12], respectively. 

It was found that f(u) is maxwellian over the entire range of current 
densities measured in [2]. The computed electrical conductivity is 
plotted as a function of current density in figure 1 along with Kerrebrock 
and Hoffman’s data. The lowest calculated value of N e was 3X10^ ra ~3, as 
noted on the figure, whereas Kerrebrock and Hoffman estimated that 
electron- electron collisions could not maintain a maxwellian f(u) if N e 
were below 3X10^-^ Neither the present theory nor Kerrebrock 1 s theory 

shows the dip in conductivity observed at 10 3 amp-nT 2 . 


Cesium atomic model 


For investigation of the effect of inelastic collisions, a five-level 
model cesium atom, which contained three discrete excited states, was 
chosen. These states (6P, P 1 ; 5D, D'; 7S) were treated as singlets, 
neglecting the slight splitting of the 6P and 5D states due to spin-orbit 
coupling. All the remaining excited states of cesium were collapsed into 
a lumped state whose degeneracy and binding energy were varied to explore 
the sensitivity of the calculation to this atomic model. Various criteria 
were employed to estimate the maximum stable electronic orbit radius. The 
features of the atomic model are summarized in Table I. 


Inelastic cross sections for electron-cesium atom collisions 


(a) Excitation-deexcitation 

Experimental [13,14] and theoretical [15,16] curves for the first 
(6S 6P) and total excitation cross sections of cesium are shown in 

figure 2. In computations the experimental [ 13] curve was used above 
2 eV, with a linear segment of slope 71 P? - eV"l [14] for 
1.42 < u < 2 eV. 

The conventional Gryzinski cross sections [ 16 ] were used for the 
remaining optically allowed transitions, and the Gryzinski exchange formula 
[17] was used for optically forbidden transitions. The classical Bohr- 
Thomson cross sections [9] provided a check on the accuracy of the compu- 
tations, since their use enables analytical evaluation of the collision 
terms for a maxwellian f(u). 

( b ) Ionization-capture 

The Gryzinski ionization cross sections [16] were employed for transi- 
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tions from all bound electronic states to the continuum. The maximum cross 
section predicted for ionization from the ground (6S) state is 20 percent 
lower than the best experimental value [18]. 


Radiative transitions 


(a) Deexcitation 

The radiative transition probabilities listed in Table I were taken 
from [ 19] . The optical thickness of the plasma can be varied simply by 
altering these values, thus simulating different plasma dimensions and 
optical absorption coefficients. 

(b) Radiative capture 

At low R e and T e , radiative capture becomes important. Cross 
sections for this process were taken from [20]. 


Solution of rate equations for maxwellian distribution 


Figures 3 and 4 show the results of solving the steady-state rate 
equations (7) for a maxwellian f(u). In figure 3 the degree of ioniza- 
tion of the cesium is plotted as a function of the initial cesium number 
density Nq s for kT e /e = 0.2 eV. The results are qualitatively similar 
to those of [4]. For an optically thin plasma, N e is determined at very 
low n£? s by the balance between radiative capture and collisional ioniza- 
tion from the ground state. With increasing N^ s , three-body capture 
becomes important, and in the optically thin case, the curve rises to a 
maximum as ionization from excited states becomes significant, thereafter 
approaching the Saha curve. For a plasma optically thick to resonance 
radiation, ionization from excited states is important at lower Nq s 
because the populations of these states are higher. Figure 3 also demon- 
strates the relative insensitivity of the results to the degeneracy and 
binding energy of the lumped state. 


Figure 4 shows that the Saha equation is a good approximation when 
T e > 3000° K, regardless of the optical thickness of the plasma. Below 
this temperature, the optical thickness of the plasma to resonance radia- 
tion is the controlling factor in determining the magnitude of the in- 
elastic term 1(f). When resonance radiation completely escapes, the 
lack of detailed balancing of collisional transitions between the 6S and 6P 
states results in a large (negative) contribution to 1(f), which is 
comparable in magnitude to the electron- electron collision term,. For a 
plasma optically thick to resonance radiation, collisional transitions 
between these states are more nearly in balance, and 1(f) is an order of 
magnitude smaller and due mainly to the 6P -» 5D and 5D -» 7S transitions. 


Distribution function and electrical conductivity - inelastic case 

Preliminary results obtained by iterative solution of equations ( 2 ) 
and (7) are shown in figures 5 and 6 for an atmospheric-pressure argon- 
cesium plasma optically thick to resonance radiation, the most realistic 
case for an MHD generator [21] . These calculations were performed for a 
seed fraction (n^/N^) of 0.0015, gas temperature of 1500° K, and no 
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magnetic field, as in the experiments of [2]. The conductivity - current 
curve of figure 5 displays the same general features as the experimental 
argon-potassium results. Of particular interest is the fact that the 
theoretical curve drops off more rapidly with decreasing current than in 
the elastic case. The reason is that the states are out of equilibrium at 
low currents and electron temperatures, and there is an appreciable in- 
elastic energy loss from the electrons, as well as a sizable reduction of 
N e below the Saha value. At the higher values of j, the two-temperature 
theory is a good approximation, as shown by the distribution function plot 
of figure 6 where the calculated f(u) at T e = 3000° K is indistinguish- 
able from a maxwellian. At 2000° K, however, f(u) departs from maxwellian 
above the first excitation potential. Since convergence difficulties were 
encountered for T e < 2300° K, the results in this range are uncertain, 
and their questionable nature is indicated by the dashed portion of the 
theoretical curve in figure 5. In the optically thin case, convergence 
difficulties were experienced at even higher values of T e . 


Conclusions 


(a) Elastic collisions cannot alone account for anomalies in the 
electrical conductivity observed at low current densities. 

(b) The Saha equation is a good approximation when either the cesium 
number density is above 1C>24 m -3 or -the electron temperature above 3000° K, 
regardless of the optical thickness of the plasma. At lower densities or 
temperatures, the optical thickness to resonance radiation controls the 
electron number density and the magnitude of the inelastic collision terms. 
Fortunately, the electron number density proves to be: relatively insensi- 
tive to the assumptions made concerning the lumped state. 

(c) Preliminary results of a self-consistent calculation of the elec- 
tron number density and distribution function show trends that agree with 
available experimental data, but more work is required to obtain reliable 
results at low electron temperatures, where inelastic effects dominate. 
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APPENDIX A - SYMBOLS 5 


B 

1 

E* 

E L 

^L,K 

e 

f(u) 


h 

1(f) 

J 


K 


L-»K 


„cap 

K L 


probability (sec - -*-) of 
spontaneous radiative 
transition L-»K 


K L 


ion 




A 


Ir>K 


K<L 

magnetic field 

electric field 

= E + U X B, effective 
electric field 

binding energy of level L, 
eV 

S E L " \ 

electron charge 

energy distribution func- 
tion (cf. footnote l) 

degeneracy 

Planck constant 

inelastic gain 

current density 

= (Se/mg) 1 / 2 /*" U f(u) 

J 0 

X Q- [r> ^(u)du, rate coef- 
ficient (m 5 - sec - -*-) 
for collision-induced 
transition L-»K 

n oo 

uf(u)u'f(u') 

X a£ a P(u' ,u)du'du, effec- 
tive two-body rate coef- 
ficient (m 5 - sec"l) for 
capture into level L 


k 

m 

N 


(u) 
ee v ' 


<o (u) 


Q£ on (u) 

T 

— * 

U 

u 

v 

W R 


5 Mks units are used for all quantities 


= (Se/mJ 1 / 2 f" uf(u ) 

e l 

X Q^ on (u)du, rate coef- 
ficient for ionization 
from level L 

$ Y - 4°" 

K 

Boltzmann constant 

particle mass 

number density 

= 4Q*.(u) 

= 4jt(e/4rte 0 )^u - ^ln A 

cross section for momentum 
transfer between elec- 
trons and species j 
(j = a,e,i) 

total cross section for 
collision- induced tran- 
sition L-*K 


-f 0 


u-E L 


a i° n ( u " jU )du", 


total cross section for 
ionization from level L 

temperature 

gas velocity 

electron energy (m v E /Se), 
eV e 


electron speed 

radiative power loss per 
unit volume 
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3l 

rate coefficient 

. 7 *i i 

Subscripts and superscripts 


(m°-sec -J -) for radi- 




ative capture into 

Ar 

argon 


level L 





a 

atom 

e o 

permittivity of vacuum 



A 

= 12rt(e 0 kT e ) 3 / 2 e' 3 N" 1 / 2 , 

Cs 

cesium 


ratio of Debye radius 

cap 

capture 


to impact parameter 




for 90° scattering 

r n 

dex 

deexcitation 


A(u) 


^ '| N h Q eh^ u ^ 


-1 


'/free 


[_h=a, i 
path for elastically- 
scattered electron 


ex 

g 


electron 

excitation 

gas 


a 


electrical conductivity H Hall 


a£ a P(u',u) probability (m 4 -sec) for 
capture of u 1 - 
electron into level L, 
with u-electron gain- 
ing excess energy 
u' + E l 

a^ on (u",u' ) differential cross sec- 
tion (m^-eV^) for 
producing a u" 
secondary electron by 
impact of a u' pri- 
mary electron 


h heavy particle (atom 

or ion) 

i ion 

ion ionization 

K potassium 

K,L atomic level indices 

r-cap radiative capture 

s seed 


T 

CD 


= A(u)/v, collision time 

= eB/m g , electron cyclo- 
tron frequency 


initial number density 
perpendicular 
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TABLE I. - PARAMETERS OF CESIUM ATOMIC MODEL 


State 

Level 

L 

Binding 

energy, 

e L> 

eV 

Degeneracy, 

SL 

Radiative transition probability, 

a k-+l> 107 sec_1 

K = 1 

K = 2 

wmm 

HKI 

K = 4 

K = 5 

6S 

1 

3.89 

2 


6.10 c 

d 

d 

0.306° 

6P,P* 

2 

2.47 

6 



0.178 

1.76 

.564 

5D,D' 

3 

2.09 

10 




d 

.378 

7S 

4 

1.59 

2 





.849 

Lumped 

5 

Varied 8, 

Varied* 5 





— 


a E 5 = 0.4, 0.6, 0.8 eV. 
b g 5 = 50, 100, 150. 

c For plasmas optically thick to resonance radiation, these are equated 
to zero. 

a Forbidden transition (A l £ +1). 





























E-3415 


Electron temperature, 



Figure 5. - Electrical conductivity as function of current density. Atmospheric pres- 
sure; gas temperature, 1500° K; seed fraction, 0.0015. 



Electron energy, u, eV 

Figure 6. - Normalized distribution function. Argon seeded with 
cesium; atmospheric pressure; gas temperature, 1500° K ; seed 
fraction, 0.0015; optically thick to resonance radiation. 


